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Abstract 

Synchronization of chaotic system may occur only when the largest conditional 
Lyapunov exponent of the driven system is negative. The synchronization with 
positive conditional Lyapunov reported in a recent paper (Phys. Rev. E, 56, 
2272 (1997)) is a combined result of the contracting region of the system and 
the finite precision in computer simulations. 
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Sensitivity to initial conditions is a generic feature of chaotic dynamical systems. Two 
chaotic orbits, starting from slightly different initial points in the state space, separate ex- 
ponentially with time, and become totally uncorrelated. As a result, independent identical 
chaotic systems cannot synchronize with each other. The sensitivity is quantitatively de- 
scribed by positive Lyapunov exponent (s) in the Lyapunov exponent spectrum of the chaotic 
system. 

However, chaotic systems linked by common signal can synchronize with each other. 
Several cases could be distinguished. In the first case, a replica subsystem driven by chaotic 
signals of the chaotic system can synchronize identically with the drive system[l-5], if the 
largest conditional Lyapunov is negative. This is referred to as identical synchronization. 

Secondly, a driven system, which is not a replica of the drive system, however, may 
not achieve identical synchronization, but generalized synchronization^-^], if the largest 
conditional Lyapunov exponent is negative. Two identical systems, driven by the same signal, 
thus may come to the same final state due to the negative largest conditional Lyapunov 
exponent. 

Lyapunov exponents are also employed to characterize behavior of random dynamical 
systems[9]: the system is chaotic (non-chaotic) when the largest Lyapunov exponent is pos- 
itive (negative). The sensitivity of a chaotic system may also be suppressed by noise, and 
identical chaotic systems subjected to common noise can synchronize with each other. Mar- 
itan and Banavar[10] studied the behavior of the noise-driven logistic maps and reported 
synchronization phenomenon. It turned out that the observed synchronization was an out- 
come of finite precision in numerical simulations [11, 12], while the Lyapunov exponent of the 
noisy logistic map is positivefll]. 

Very recently, Shuai et ai[13] claimed that synchronization can be achieved with positive 
conditional Lyapunov exponents. In a one-way coupled map lattice, they observed, through 
computer simulations, synchronization of spatiotemporal chaos with many positive compo- 
nents in the conditional Lyapunov exponent spectrum. Based on these results, they drew 
the conclusion that the conditional Lyapunov exponents cannot be used as a criterion for 
synchronous chaotic systems. 

Whether such a claim is true is of great importance for our understanding of synchro- 
nization. In this paper, we reexamined such synchronization phenomenon, revealing that it 
is yet another example of round-off induced phenomenon. 

In[13], Shuai et al studied a driven one-way coupled map lattice 



Vi (t + 1) = (1 - e)f( yi (t)) + ef(y i+1 (t)) (i = 1, • • • , N), 

VN+lif) = X (t), 

where xo(t) is a hyperchaotic signal from a one-way coupled ring lattice 



(1) 
(2) 



Xi (t + 1) = (1 - e)f( Xi (t)) + ef{x i+1 (t)) (i — 1, ■ ■ • , N), 

x N+1 (t) = x (t), 




(3) 
(4) 
(5) 



The chaotic map is the well-known logistic map f(x) = 4x(l — x) and eo = 0.01 
As in [13], the conditional Lyapunov exponents of the driven system are 



A i = ln(l-c)+ Km ^\f( yi (t))\ 




(6) 
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Let us study the simplest case of N = 1. The conditional Lyapunov exponent as a 
function of e is shown in Fig. 1. To detect the behavior of synchronization, 100 performances 
with random initial conditions are carried out for each e. Synchronization occurs when y 1 
and Xi become numerically identical for the finite precision in simulations (double precision). 
P = M/100, where M is the number of simulations in which synchronization occurs within 
5 x 10 7 iterations, is estimated as a function of e, as shown in Fig. 1. It can be detected 
that synchronization with positive Lyapunov exponent (SP) occurs in several regions. We 
will take e = 0.200 and e = 0.335 as examples, as was pointed out in [13]. 

Is SP a true physical phenomenon or an artifact of finite precision in computer simula- 
tions? In the following, different precision formats (single, double and quadruple precision) 
are employed in the simulations. Firstly, the difference e(t) = \yi(t) — xi(t)\ preceding the 
synchronized state is examined for simulations with different precisions but the same ran- 
dom initial conditions. The results for e = 0.200(Ai = 0.105) and e = 0.335(Ai = 0.025) are 
shown in Fig. 2(a) and (b), respectively. For comparison, an example of synchronization 
with negative conditional Lyapunov exponent (SN) at e = 0.520(Ai = —0.041) is illustrated 
in Fig. 2(c). Note the different scales for the different precisions used. e(t) displays an 
intermittent behavior before reaching SP. SP occurs somewhat abruptly, when e(t) drops 
lower than the precision of the computer. The time T needed for SP to occur is much longer 
for quadrupole precision than that for single and double precisions. As for SN, e(t) continues 
its trend of decrease when higher precision is employed in the simulation. It is plausible to 
imagine that for SP, the intermittence of e(t) will continue indefinitely for infinite precision, 
while for SN, e(t) will approach to 0. So, the physical process of SP is an intermittence, 
with e(t) becoming very small and enlarging to the size of the chaotic attractor alternately. 

The difference e(t) is actually not zero even beyond the precision of the computer. In the 
following simulation with quadrupole precision, when the states of the systems are numeri- 
cally identical (SP), a perturbation £ e (— 10~ 30 , 10 -30 ) is added to the drive signal x (t) of 
system y at the next iteration, under the constrain < xo(t) + £ < 1. Such a tiny perturba- 
tion can totally destroy the synchronization behavior when Ai > 0, as seen from Fig. 3(a) 
and (b) for the results of e = 0.200 and e = 0.335 respectively, because the tiny difference 
can be amplified to the order of 10° due to the positive conditional Lyapunov. While for 
SN, e(t) continues to decrease after the impulsive perturbations, and the level of difference 
is the order of 10 -30 . Such a dramatic difference between the behavior of SP and SN shows 
that, negative conditional Lyapunov is necessary condition for physical synchronization. 

To demonstrate further that synchronization cannot be achieved physically with a pos- 
itive conditional Lyapunov exponent, the average synchronization time T a is evaluated for 
different precisions with 100 random initial conditions. The results for SP (e = 0.200 and 
e = 0.335) are displaced with a linear-log plot in Fig. 4(a). The three points lies almost 
on a straight line, meaning that T a ~ exp(AL), where L is the number of significant digits 
of the finite precision. An exponential increase of T a with L proves that synchronization 
can never occur with infinite precision. The behavior of SN is greatly different, where T a 
follows a linear dependence on L, T a ~ BL, as seen from the result of e = 0.52 displaced 
with a linear-linear plot in Fig. 4(b). The reason is that, approximately, e(t) decreases 
exponentially with time, so that 10 _L ~ exp(AiT a ), resulting in B — — lnl0/Ai. B = 56.2 
at e = 0.520 is in good agreement with the slope 58.5 of the solid line in Fig. 4(b). 

Why is then that SP can be observed in numerical simulations even within thousands of 
iterations? The origin is that there are contracting regions in a map f,C— {x,y\\f(x) — 
f(y)\ < !}• Two orbits in a contracting region come closer to each other at the next step. 
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For the system studied above, the contracting region is 1 — jfjz^j < x i + Vi < 1 + 4(i- £ ) • 
The strip near xi + t/i = 1 has the strongest contracting rate. The distribution of x\ + yi 
is calculated with 10 7 iterations to examine the relationship between SP and C. As seen 
from the results of e = 0.200 and e = 0.335 in Fig. 5(a), the distribution for e = 0.200 has 
very high peaks in the contracting region, while a lower peak for e = 0.335. Such greater 
frequency of access to the contracting region at e = 0.200 makes e{t) drops to a much small 
value more frequently than at e = 0.335 (see Fig. 3), which accounts for the result that 
SP is observed at e = 0.2 with much fewer iterations than that at e = 0.335. However, the 
evolution of the difference e(t) is a combined result of local stability and instability. The 
finite-time Lyapunov exponent [14] 

i m 

X^ = -^\f'( yi (t))\ (7) 

measures the average expansion or contraction rate in m steps. The distributions of X^ (m = 
70) for e = 0.200 and e = 0.335 are illustrated in Fig. 5(b). A pronounced tail to —1.0 at 
e = 0.200 means that the difference shrinks by a factor of e~ 70 = 4 x 10~ 31 in some successive 
70 iterations. The negative tails thus plays an important role in observation of SP. It also 
explains the fact that SP is easier to occur at e = 0.200 with a larger positive conditional 
Lyapunov exponent (A = 0.105) than at e = 0.335 (A = 0.025). The distributions of A^ 
also reflects the true dynamics of SP: the difference e(t) can be very small in a period of time, 
and it will be amplified in some other period of time because temporal separation dominates, 
thus resulting in an intermittent dynamics. So the finite-time Lyapunov exponent gives a 
more convincing account for the occurrence of SP in simulations. 

Based on the above analysis of the simplest case of N = 1, the SP observed in computer 
simulations is a combination of two factors: the shift of the state of the chaotic system 
to the contracting region and the finite precision in numerical simulations. For the case 
of N > 1, SP is observed for similar reasons. The driven system is coupled in a cascade 
way x N — > Xn-i — > ••• —> %i, and synchronization can only occur for the first several 
Ni nodes if all the Ni conditional Lyapunov exponents Ajv, Ajv-i, • • • Xn-Nx+i are negative. 
Physical synchronization of all the lattices can only occur when all the conditional Lyapunov 
exponents are negative. 

In conclusion, synchronization with positive conditional Lyapunov exponents in com- 
puter simulations is a round-off induced phenomenon. The physical dynamics of SP is an 
intermittence. One can expect to observe SP easily in computer simulations in such systems 
with large contracting regions, and the couplings have the effect of shifting the state to the 
contracting regions so that the finite-time Lyapunov has a significant tail of negative values. 
A negative conditional Lyapunov exponent is a necessary condition for synchronizing chaotic 
systems. 
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Figure Captions 



Fig. 1. Conditional Lyapunov exponent Ai and synchronization ratio P as functions of e. 

Fig. 2. The time series of the difference e(t) proceeding synchronization in simulations 
with single, double and quadrupole precisions. (a)e = 0.200, SP; (b) e = 0.335, SP; 
and (c)e = 0.52, SN. 

Fig. 3. The time series of difference e(t) under impulsive perturbations between (10 -30 , 10 -30 ) 
in simulations with quadrupole precision. The initial conditions are the same as Fig. 
2. (a)e = 0.200, SP; (b) e = 0.335, SP; and (c)e = 0.52, SN. 

Fig. 4. Average synchronization time T a as a function of precision in simulations. L — 7, 16 
and 31 for single, double and quadrupole precision, respectively. The each solid line 
links the first and last point of each data set. (a) Linear-log plots for SP. (b)A linear- 
linear plot for SN. 

Fig. 5. (a) Normalized histograms of x\ + y\. (b) Normalized histograms of finite-time 
Lyapunov exponent. 
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